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TECHNICAL MEMORANDUM 


ANALYTICAL ROUND ROBIN FOR ELASTIC-PLASTIC ANALYSIS OF SURFACE 
CRACKED PLATES, PHASE II RESULTS 


1. INTRODUCTION 


The ASTM Committee on Fatigue and Fracture (E08) Task Group on Fracture Toughness 
of Surface Cracks (E08.07.03) has developed a new material test standard for the assessment of 
surface crack toughness in the linear-elastic and the elastic-plastic regime, E2899.! While evalu- 
ation of the linear-elastic stress intensity factor, K, is well-accepted for surface crack geometries 
through the Newman-Raju equations” or other existing tabulations,‘ the evaluation of elastic- 
plastic J-integral values requires the user to perform an analysis of the experiment using methods 
that cannot be easily contained or conveyed within a test standard. The feasibility of allowing this 
latitude in assessment methodology in a testing standard requires validation through an interlabo- 
ratory study (ILS). This Technical Memorandum (TM) provides an overview of the problem state- 
ment and results of the second phase of the study. 


1.1 Motivation for the Interlaboratory Study 


Existing mechanical testing standards for fracture mechanics applications, such as ASTM 
E1820,° include within the document all the necessary equations and supporting data to assess the 
experiment directly, producing a ‘standard’ result. In the case of a surface crack in a flat plate, no 
openly available method exists for accurately evaluating the J-integral under elastic-plastic condi- 
tions along the full perimeter of a surface crack that can feasibly be reduced to equation form for 
inclusion into a test standard. However, the ability to assess this problem using methods of analysis 
such as finite elements has progressed considerably in recent years and may now be a suitable substi- 
tute for assessment of experimental results in lieu of equations defined within the standard. Codes 
defining standard practice for structural evaluation of defects, such as the API Recommended Prac- 
tice for Fitness for Service, allow for such external analysis methods. This study is designed to eval- 
uate the variability in the use of an external analysis method for experimental data assessment. 


1.2 Scope and Objectives for the Interlaboratory Study 


This ILS is planned for development in phases based on need and observations made as 
the study progresses. Additional phases are likely to investigate different crack shapes, bending 
as opposed to tension loading, and different materials. This work was conducted under the aus- 
pices of ASTM ILS 732. The first phase of the ILS is documented in NASA/TM—2012-217456,° 
and approached the ILS concept differently than most experimentally based studies. The analysis 


methodology was not specified to participants and key experimental results were withheld. This 
made phase I of the study ‘blind’ in the sense that the participants did not have the experimental 
force and displacement information available to them to validate their model, or to potentially alter 
their model to converge on the experimental result, and they were free to approach the analysis 
using any method of choice. This approach allowed this first phase of the ILS to serve as a cur- 
rent measure of the state of the art for elastic-plastic fracture mechanics analysis. An open study 
of this nature regarding elastic-plastic analysis methods had not been performed for some time.” 
Examination of the phase I results, as detailed in reference 6, isolated the sources of variability in 
the finite element analysis (FEA) and identified the common analytical practices. These findings led 
to the contents of Annex A6 in E2899, “Methodology for performing elastic-plastic finite element 
analysis and comparison with test record.” 


The objective of phase II is to enhance the findings from phase I of this activity by continu- 
ing to collect independent analysis of surface cracks in tension by requesting participants perform 
an FEA of a new experiment while adhering to the methodology of the E2899 Annex A6 as closely 
as possible. The expectation being that, if each participant’s analysis follows the best practices 
defined by Annex A6, the variability in results would be reduced relative to the (already favorable) 
observations in phase I. This second phase also differs significantly from the first phase in that the 
full experimental test records are provided to the participants. In addition to the FEA, the study 
participants were asked to interpret the test results according to E2899. This involves evaluating 
crack-front constraint, determining a critical initiation angle (@;) along the crack front, determining 
the deformation regime of the test (linear-elastic fracture mechanics (LEFM), elastic-plastic fracture 
mechanics (EPFM), or field collapse), and determining the critical J-integral value. The intent was 
that this process would familiarize participants with the standard and provide the task group with 
helpful feedback regarding the clarity of the language in the written requirements of the standard. 


In addition to the FEA of the test, the task group is interested in evaluating a potential new 
means of analysis of surface crack experiments suitable to the standard. Allen?:!° has developed an 
interpolation tool to evaluate surface cracks in tension across a wide range of semi-elliptical crack 
aspect ratios and material behavior. The interpolation tool works with a broad database of 600 
solutions of surface cracks under a deformation state extending to the limits of the elastic-plastic 
regime. The software that interpolates this space, Tool for Analysis of Surface Cracks (TASC), 
is available free of charge.!! For participants willing and able to perform FEA but also interested 
in TASC, the authors requested that the FEA work be completed prior to assessing the test with 
TASC. For those interested in evaluating the test but unable to participate in the FEA work, study 
participants were also allowed to evaluate the problem using only TASC. 


The experimental test for this second phase of the ILS was planned as a ductile tearing 
test on steel with moderate yield strength providing a result well into the EPFM regime, beyond 
the limits of LEFM validity. The choice of readily available 4142 steel proved inappropriate when 
unstable cleavage failure occurred—it had not occurred to us that the transition temperature for 
4142 may be as high as 100 °C. Our intent was to produce a fracture surface with an identifiable 
location of maximum stable ductile tearing along the crack front so participants could identify 
the critical initiation angle from an image of the fracture surface. This would provide insight into 
the repeatability of this required image measurement. In this case, though the test did reveal some 


stable tearing, the fracture surface did not provide obvious evidence in the ductile tearing of the ini- 
tiation angle for test evaluation. Instead, participants were requested to work through the process 
described in E2899 Annex AS to identify the initiation angle by using the product of the J-integral 
and constraint (7-stress). For this ILS activity, the authors chose to proceed with an assessment of 
this test because it fully suits the purpose, in spite of section 1.8 of E2899 which declares cleavage 
results being outside of the scope of the standard. 


2. PROBLEM STATEMENT 


The participants were asked to evaluate a surface crack tension test according to the version 
of the standard published at that time, E2899-13. The basic geometry and definition of parameters 
for a generic surface crack specimen are shown in figure 1. The complete problem statement as it 
was provided to each participant is included as appendix A. A compressed archive file was provided 
to each participant containing the following data: instructions, material identification, and uniaxial 
tension stress and strain data, force versus crack mouth opening displacement (CMOD) record for 
the surface crack tension test, and the surface crack fracture surface photo with basic crack size 
measurements. These data represent the fundamental information provided by a mechanical test 
laboratory prior to evaluating the test result. The choice of FEA constitutive models and treatment 
of the stress-strain data was left to the discretion of the participants. 


From their evaluation of the test per E2899, each participant was requested to provide sev- 
eral results including the fracture toughness test result in K or J at the critical initiation angle (@,), 
crack front constraint quantification in terms of 7-stress, and interpretation of selected validity 
criteria. Spreadsheets for recording analysis results were provided to enhance consistency in report- 
ing and to clarify what data were being requested. See appendix A for details. 
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Figure 1. Surface crack test specimen configuration. 


3. PARTICIPANT OVERVIEW 


This phase of the study had 10 participants covering a range of industries, academia, and 
nationalities. There was also intent to include participants with a broad range of expertise in the 
study. The user community for the test standard is potentially broad; therefore, our intent was not 
to have only ‘experts’ in the field participating. A study with only expert participants would not be 
properly representative, leading to overly optimistic results and a lack of proper feedback on the 
test standard. 


A list of the participants is included in table 1. All participants have been assigned a random 


designation from Lab-1 to Lab-10, with Lab-1 representing the authors’ results. All results remain 
anonymous, though each participant was informed of their lab identification. 


Table 1. ILS participants and affiliation. 


Name Affiliation 


Phillip Allen MSFC Materials Laboratory 

Steven Altstadt Stress Engineering Services 

Jason Bely Alcoa 

Enrico Lucon National Institute of Standards and Technology 


Francisco Martin | Purdue University 

Dawn Phillips MSFC Vehicle Structures 

Ryan Sherman Purdue University 

Greg Thorwald Quest Integrity Group 

Igor Varfolomeev | Fraunhofer-Institut fuer Werkstoffmechanik IWM, Germany 
Michael Windisch | MT Aerospace AG, Germany 


4. RESULTS AND DISCUSSION 


The phase II experiment consisted of a 4142 steel surface crack tension specimen of the 
dimensions shown in appendix A. A photograph of the specimen during test (fig. 2) shows the ring 
gauge used for CMOD measurement. The specimen was precracked in tension, instrumented, and 
then loaded monotonically while monitoring the force and CMOD values. The test intent was to 
load the specimen to an estimated critical CMOD value for the initiation of ductile tearing and 
then unload and fatigue marker-band the specimen to highlight any region of crack tearing. But, 
as mentioned earlier, the specimen failed in an unstable fashion due to cleavage before the esti- 
mated critical CMOD was reached. After consideration, the decision was to not repeat the test 
but to proceed with this result for the ILS. 


Figure 2. ILS specimen configured for testing. 


The fracture surface for the ILS test (fig. 3) shows a largely uniform ductile crack extension 
(¢) of approximately 0.25 mm. Note that the cleavage initiation site can be identified, but this was 
not mentioned to the ILS participants. The minimum precrack extension was 1.30 mm and the 
electrical discharge machined initial notch height, NV, was 0.43 mm. The CMOD value at specimen 
failure was 0.196 mm corresponding to a failure force of 603.5 kN. This information, combined 
with the force-CMOD response, constitute the experimental test result. 


Figure 3. ILS specimen fracture surface. 


4.1 Tying Analysis to Experiment 


For a test result that falls outside the LEFM regime, the only physical tie between the test 
analysis and experiment is the force versus CMOD response. Based on this metric, E2899 requires 
agreement between the analysis and the experiment in two ways:—the slope of the force-CMOD 
elastic response and the final analytical force-CMOD prediction. 


The force-CMOD response from the experiment is shown in figure 4. The nonlinearity 
in the test record is an indicator of elastic-plastic conditions in the specimen, and no unloading 
slope is plotted since the specimen exhibited unstable fracture. The force-CMOD predictions from 
the 10 ILS labs are plotted along with the experimental record in figure 5. Note that since some 
labs performed the analysis using both FEA and TASC, there is a total of 12 predictions, and the 
TASC results are designated by a ‘—7” suffix to the lab designation in the legend. The analytical 
predictions of the force-CMOD response show good overall agreement. Specific evaluations of the 
force-CMOD predictions in the linear region and at the failure CMOD of the test are discussed 
separately. 
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Figure 4. Experimental force versus CMOD response. 
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Figure 5. Force versus CMOD response for all ILS participants. 


E2899 does not specify a procedure to determine the slope of the linear portion of the 
experimental force-CMOD response, so the slope of the linear portion of the test response was 
estimated using two methods: (1) Linear regression of the force-CMOD data falling within a fixed 
range of 20%-50% of the maximum force, and (2) the slope determination by analysis of residu- 
als (SDAR) method as developed by Graham and Adler.!? The elastic slope of each analysis 
result was determined using the first two load steps in the analysis where the material response 
is predominantly linear. The linear regression slope over the fixed range was 8,112 kN/mm and the 
SDAR slope was 8,291 kN/mm, for a difference of 2.2%. This reflects the long-standing concern 
regarding the challenge of accurately determining the slope of experimental records. The percent 
difference between the determined elastic slopes of the experiment and each analysis is presented 
in table 2. For the fixed range slope, all analytical predictions were within +5% except for Lab-2-T 
and Lab-4. Using the SDAR method, all analytical predictions were within +5% except for Lab-4 
and Lab-9-T. For reference, Annex A6.3 of the E2899-13 version used by the ILS participants 
specified that “the most linear portion of the force-CMOD test record and elastic compliance of 
the finite element analysis result shall match within 2.5%.” Using this metric, only about half of 
the participants would have met the requirements for a valid analysis result, but, as will be shown 
later, the J-integral results between the labs are still in close agreement. The elastic compliance 
check in E2899 is intended to serve as a broad check on numerous basic inputs to the analysis such 
as specimen and crack geometry, material elastic properties, and boundary conditions to screen out 
fundamental analysis mistakes, but was not intended to be so stringent as to preclude otherwise 
acceptable analysis results. After studying the results of this ILS and considering the expected accu- 
racy of slope determination from the experimental record, an elastic slope matching the require- 
ment of +5% was a more reasonable metric to meet the intent of the elastic compliance check. This 
change was implemented in the E2899-15 release of the test standard. 


Table 2. Summary of elastic compliance evaluation results. 


Experiment Elastic Slope Experiment Elastic Slope 
Determined Using Linear Determined Using SDAR 
Fit to 20%-50% of Max Data Range |) Graham-Adler Fitting Algorithm 


Lab Elastic Slope % Diff. Elastic Slope % Diff. 


Lab-1 0.06 2.23 
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The second metric needed to correlate the analysis to the experimental result involves the 
values of force and CMOD recorded at the initiation of crack tearing, or, as in the case of this test 
result, at the point of unstable fracture. Since it is unlikely the test and analysis records will match 
exactly, an analyst must choose to evaluate the test result based on the analysis matching either the 
test initiation force, P,, or the test initiation CMOD, CMOD,. Though force is an intuitive descrip- 
tor of the test failure point, for elastic-plastic analysis, the CMOD is a more reliable predictor of 
the J-integral, because J is nearly a linear function of CMOD in the plastic regime.® !9!3.14 E2899 
specifies that the user determine the analysis reaction force corresponding to an analytical CMOD 
equal to CMOD, from the experiment. For an analysis result to be considered valid, the analysis 
reaction force corresponding to CMOD, shall be within +5% of the experimental P,. The green star 
symbol and error bar in figure 5 illustrates the test evaluation point and the +5% force error band 
for this experiment. All of the ILS analytical force results corresponding to CMOD, easily fall 
within the +5% force error band except for Lab-9-T. The Lab-9-T force value exceeds P; by 5.25%, 
but, since the results are just outside the chosen bounds, the Lab-9-T results are still included in the 
ILS. 


4.2 Fracture Mechanics Analysis Results 


To begin an analysis of a surface crack test, section 9.2 of E2899 requires users to determine 
the deformation regime of the test: linear-elastic, elastic-plastic, or field collapse regime. Determin- 
ing the deformation regime is fundamental because each regime requires different treatment of 
the data for valid results calculation and reporting. Per E2899, for a test to be in the linear-elastic 


regime, the net section stress calculated using P; must be less than 90% of the yield stress, o,,,, and 
Toa> "pb 2Cr Ix [Fys); where ; 
(Ko) (I-v") 
Fee (1) 
E 
and 
E 
Cc ‘Ko Oo. ’ (2) 
ys 
where 
Kg = linear elastic stress intensity factor calculated using P, at the initiation angle, @; 


v = Poisson’s ratio 
loa and = characteristic lengths of the specimen test section at @,. 


Figure 6 illustrates the definition of these characteristic lengths. Considering a 2D slice through the 
specimen along a path defined by "ga and rgb then ae is analogous to a crack length measurement 
and ry, relates to a remaining ligament measurement. For the ILS phase 2 test, all of the partici- 
pants correctly determined that both linear-elastic validity criteria were violated and, therefore, 

an elastic-plastic fracture mechanics analysis of the test must be performed. 
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Figure 6. Cross section through the crack plane illustrating the characteristic 
lengths ry, and rg), 


It is interesting to note that a test such as this, where @; is not clearly identifiable in the 
fracture surface, leaves the user in a bit of a conundrum. The user needs to know @; to be able to 
calculate Ky, rg,, and rg,, but as described in Annex A5.2 of E2899, an elastic-plastic analysis is 
required to determine @,. In practice, a user can quickly perform the linear-elastic net section stress 
check, and if it is violated, move on to the more complicated elastic-plastic analysis step. 


The as-reported J-integral versus @ results at CMOD, from all ILS participants is shown in 
figure 7, where @ follows the common convention as defined in figure 1, with @=0 at the free sur- 
face and @= 90 at the centerline crack depth location. Recall that, for this test, CMOD, occurred 
at the point of unstable fracture and is therefore defined as the maximum CMOD measured prior 
to failure. In the initial review of the analysis results, the clear outlier was that the as-reported 
Lab-8 J-integral results were approximately one-half of the other participant’s results. The Lab-8 
reporting error is a common mistake in finite element domain integral calculations wherein the 
user forgets to double the finite element domain integral output due to the symmetry plane on the 
domain boundary. Once confirmed as the error, the Lab-8 J-integral results are included, thus they 
have been properly doubled for all reported values in this TM. This instance illustrates the value of 
an ILS exercise to highlight common mistakes so they can be addressed and prevented. Even after 
correcting the Lab-8 results by doubling them, it was observed that the Lab-8 results along with 
the Lab-5 results do not agree with the rest of the participant’s results. With some investigation, it 
was discovered that the J values reported by Lab-8 and Lab-5 correspond to the analysis step where 
the FEA reaction force equals P; instead of reporting J values at the FEA step corresponding to 
an analytical CMOD equal to CMOD, as required in E2899. By anchoring their analysis result on 
force rather than CMOD, Lab-5 and Lab-8 reported their results at CMOD values smaller and 
larger than CMOD,, respectively, resulting in J-integral distributions below and above the rest of 
the result family. With the data provided by Lab-5 and Lab-8, it was not possible for us to deter- 
mine J-integral values corresponding to CMOD, for each @ location from their results; therefore, a 
final consistent J-integral versus @ result at CMOD, without Lab-5 and Lab-8 values are shown in 
figure 8. The overall shape of the J-integral distribution and range of J-integral values are a close 
match for most of the participants. Lab-4’s result has a different behavior at ¢=0, likely due to 
numerical free surface effects, and underestimates the increase in J values in the 15°< @< 45° range. 
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Figure 7. J-integral versus @ at the critical CMOD as reported by all participants. 
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Figure 8. J-integral versus @ at the critical CMOD. 


The next step in the fracture mechanics assessment of the test is the estimation of the criti- 
cal initiation angle. As mentioned before, since @; could not be visually determined from ductile 
tearing evidence on the fracture surface, E2988 directed ILS participants to Annex AS for the @; 
estimation methodology. Technical background information on the development and use of the 
E2988 Annex A5 methodology for @; estimation is given in appendix B. The following provides 
a brief summary. If @; cannot be determined from the fracture surface, Annex A5 directs the user 
to estimate @; by finding the value of @ which maximizes /(@): 


p= 40{ E,,| for TO) <9 


Jy ys o ys 
G@je1 7 Os J for £59, 3) 
where 
J() = J-integral as a function of @ 
J = peak J-integral value for all 


p 
T(o)/ Oy, = T-stress as a function of @ normalized by the yield stress 
(tables are given in E2899 Annex A2). 


Numerous previous studies have shown that, in most cases, initiation under uniform tension stress 
does not occur at the maximum value of the J-integral along the crack perimeter location, imply- 
ing that the crack front constraint conditions are influencing the angle of initiation.!>-!8 Therefore, 
equation (3) estimates the crack front location where the driving force and constraint combine to 
provide the highest likelihood for crack extension, based on the value of @ that maximizes their 
product. Figure 9 reproduces figure A5.2 from E2899 illustrating the methodology for the @; 
estimation using equation (3) (labeled as equation (A5.2) in the figure). 
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Figure 9. Figure A5.2 from ASTM E2899 illustrating the critical angle evaluation. 


The calculated @; for all of the ILS participants is given in table 3a. Lab-10 had an error in 
their 7-stress calculation which resulted in an incorrect @; value of 25° as reported. Using the cor- 
rect T-stress, the corrected Lab-10 @; value is 35° which falls in line with the average of the results 
from all participants as shown in table 3b. Even including the Lab-10 results, the critical angle 
calculation results are very consistent with less than a 5° standard deviation in the result. Figure 10 
shows the ILS results plotted on a sketch of the crack geometry with the predicted values shown 
by red markers on the crack front perimeter. Due to symmetry, only the 0°< @<90° results are plot- 
ted but another set of critical angles exist and could be plotted in the left-hand quadrant as well. 
Figure 11 shows the same @; predictions plotted on the fracture surface of the ILS specimen. It is 
noteworthy that, to the trained eye, the cleavage river patterns on this fracture surface point back 
tothe predicted initiation location, which serves as another affirming validation of the methodology. 
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Table 3a. Summary of critical angle Table 3b. Summary of critical angle evaluation results 
evaluation results. with corrected 7-stress for Lab-10. 


6; (deg) g; (deg) 


Max. 42.8 

: 
Min. 25.0 in. 
Std. Dev. Std. Dev. 


Figure 10. Illustration of the critical angle values for all ILS participants plotted 
on a sketch of the crack front. 
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Figure 11. Ilustration of the critical angle values for all ILS participants plotted 
on the test specimen fracture surface. 


With the initiation angle established, the critical J-integral value (J, o) corresponding to 
CMOD, and occurring at @; can be determined. Figure 12 shows the J-integral versus CMOD 
results at @; for all participants, and figure 13 gives a closer view of the results at CMOD equal to 
CMOD,, In addition, the individual Jy values are given in table 4. Each result set is taken at the @; 
location calculated by each individual participant shown in table 3a. Since the authors did not have 
the complete model results for each ILS participant, it was not possible to evaluate the J-integral 
versus CMOD results at a single, consensus @; value. Even so, the J-integral versus CMOD results 
are in excellent agreement for all of the analyses with <2.5% standard deviation from the average 
J result. In addition, all of the J 6 results are within 5% of the average J, 6 
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Figure 12. J-integral versus CMOD. 
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Figure 13. Closer view of the range of J-integral values at the critical CMOD. 


Table 4. Summary of J-integral values at CMOD,. 


J (kJim2) at ¢; 


The final step in the elastic-plastic fracture mechanics analysis of the test is the assessment 
of the deformation regime to determine if the J, results are valid per the deformation validity 
criteria in E2899 by evaluating the specimen characteristic lengths, ry, and rg), against their respec- 
tive elastic-plastic regime limits. Both evaluations must be true to report a valid result in the elastic- 
plastic regime. First, the amount of crack-tip opening displacement must be a small fraction of the 
crack size such that rg,2 Cj(Jg/0,,) where 


G15: (4) 


Second, the remaining ligament must have a sufficient size relative to the deformation such that 
rgp2 Cip(Jg/O,s) where 


1 E 
Cy, =——— +50. 5 
Jb~ 99 Ors ©) 


All of the participants correctly reported that the J, 6 results of this test are valid elastic-plastic 
regime results. 


Though not necessary to perform the assessment, a graphical representation of the defor- 
mation regime and limits is shown in the Lab-1 TASC-generated plot in figure 14 to provide better 
insight into the deformation level relative to the E2899 EPFM limits. Figure 14 is analogous to fig- 
ure 8 in E2899-13 (given as fig. 15 in this TM) wherein the vertical axis is the inverse of the defor- 
mation level (1/C;) to have an increasing deformation plot in an increasing sense and the horizontal 
axis is the normalized crack front constraint, here T-stress/o,,,. In figure 14 the blue and red curves 
represent the deformation trajectories for rg, (crack depth) and lob (remaining ligament), respec- 
tively, at the calculated critical angle of @;= 36° while the gray curves comparatively show the same 
assessment for 6=90°. The deformation trajectories start at (0,0) and proceed upward and to the 
left as negative T-stress and J increase with increasing deformation. The green stars represent the 
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point of unstable fracture in the test at CMOD=CMOD,. The inverse of the deformation limits for 
LEFM (1/Cx), EPFM remaining ligament (1/C7,), and EPFM crack depth (1/C,,) are plotted as 
black, red, and blue horizontal lines, respectively. For an assessment to be within the validity limits, 
the assessment point for the analysis must fall on or below its respective deformation limit line. The 
assessment points for this test are clearly above the LEFM limit, but both the remaining ligament 
and crack depth assessment points are well below their respective EPFM limits, illustrating the 
validity of the Jy values per E2899. 


Tearing Load 


T-Stress/o,, 


Figure 14. Lab-1 TASC assessment of crack front deformation conditions. 
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At initiation of ductile tearing in a test sample or structure, the crack tip = 
J conditions will fall into one of the five regions A-E in the constraint/deformation |-.- 
= diagram below. Evaluate the constraint (W) and the deformation limits (C) 
ys at the onset of ductile tearing to determine the applicable region for 
assessment of crack tip conditions. 


©) ct 


Lelelelelelelelelebelelelslelelelslelelelelelelelelelelse LL Lt fe 1IC 


Cie 


Increasing Deformation ——»> 
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required to describe fields 0 one parameter required 
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Figure 15. Figure 8 from ASTM E2899 illustrating the crack front deformation 
conditions and test regime assessment. 


4.3 Comparison of Finite Element Modeling and the Tool for Analysis of Surface Cracks Results 


In addition to the primary focus of phase II of the ILS, there was high interest in evaluating 
a potential new means of analysis of surface crack experiments suitable to E2899 by using interpo- 
lated solutions and the TASC software tool. The advance of surface crack fracture testing is hin- 
dered significantly by the lack of a readily available set of solutions to correlate the applied force 
and observed CMOD in a surface crack experiment to an evaluation of the elastic-plastic J-integral 
or deformation state of a test specimen at fracture. As reflected by the necessity for this ILS, the 
only practical way to fully analyze such a test is through the use of elastic-plastic FEA. A conve- 
nient and practical set of elastic-plastic surface crack solutions could help mitigate many of these 
obstacles; however, to date, it has proven impractical to reduce the 3D elastic-plastic surface crack 
solution to a set of equations suitable for inclusion in a testing standard. 


To address this issue, the authors developed the computer program TASC with a graphical 
user interface that allows easy access to a comprehensive space of nonlinear J-integral solutions 
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for surface cracks in tension. TASC employs a methodology for interpolating between the geomet- 
ric and material property variables that allows the user to estimate the J-integral solution around 
the surface crack perimeter (@) as a function of loading condition from the linear-elastic regime 
continuously through the fully elastic-plastic regime. In addition to the J-integral solution, the 
complete force versus CMOD record is estimated to provide a direct anchor to the experimental 
result. The user of this interpolated solution space need only know the crack and plate geometry 
and the basic material flow properties to reliably evaluate the full surface crack J-integral and force 
versus CMOD solution. Thus, a solution can be obtained very rapidly by users without elastic- 
plastic fracture mechanics FEM experience. 


Four ILS participants chose to use TASC to evaluate the test result—two were experienced 
analysts that used TASC in addition to their nonlinear FEA assessment (Lab-1 and Lab-2) while 
two other participants had limited surface crack fracture test analysis experience (Lab-6 and Lab-9) 
and used only TASC. In addition, Lab-9 used TASC’s default material stress-strain curve fits to the 
provided data for the nonlinear material property inputs in the analysis without realizing that it is 
acceptable and desirable to adjust the default curve fit values to better fit the actual stress-strain 
response if needed. Even so, the TASC assessment results are clearly in family with the nonlinear 
FEA results. Comparing both the ¢; and J, values from the FEA and TASC assessments shows 
that the average results are within 2%, illustrating that the TASC assessments provide an answer 
of equivalent quality to that expected from a custom finite element assessment of the test. 
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5. PARTICIPANT COMMENTS 


The participant comments from phase II of the ILS were encouraging. In general, the 
participants were able to follow the text of E2899 and work through the assessment without issues, 
obtaining results with little variation. Some typographical errors and confusing language were 
pointed out by participants which resulted in a number of changes in the E2899-15 version of the 
standard. By far the most common technical difficulty was in matching the initial elastic slope of 
the analysis and test record to the E2889-13 imposed slope matching limit of +2.5%. As mentioned 
earlier, after studying the results of this ILS it was determined that an elastic slope matching the 
requirement of +5% was a more reasonable metric to meet the intent of the elastic compliance 
check, which is reflected in the E2899-15 release. 
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6. CONCLUSIONS 


Phase II of the ILS was fully successful with encouraging results. Ten participants with 
a wide range of experience volunteered to analyze a nonlinear surface crack test result following 
the guidance in E2899-13. The guidance on performing elastic-plastic FEA developed from phase I 
of the ILS and incorporated into E2899-13 provided sufficient instruction, allowing the partici- 
pants to independently perform custom elastic-plastic FEA of the test. The final J, test results were 
in close agreement, certainly well within expected errors in other nonlinear fracture toughness test 
methods. In addition, participants that chose to use TASC to perform the analysis achieved results 
equivalent to those performing the custom FEA. Participant comments from phase II were evalu- 
ated and incorporated into changes in the E2899-15 release, resulting in an improved test standard. 
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APPENDIX A—PHASE IT PROBLEM STATEMENT 


ASTM Committee on Fatigue and Fracture 
E08.07.03 Task Group on Surface Cracks 


Analytical Round Robin on Elastic-Plastic Analysis 
of Surface Cracks in Flat Plates — Phase II 


Introduction: The £08.07.03 task group published ASTM E2899 - Standard Test Method for Initiation 
Toughness in Surface Cracks Under Tension and Bending in October of 2013. This standard enables fracture 
toughness testing in linear-elastic and elastic-plastic regimes. In Phase | of this round robin, participants 
were asked to evaluate a surface crack in tension using finite element analysis. Those results can be found 
in NASA-TM-2012-217456, which can be retrieved from the NASA Technical Reports Server at 
http://ntrs.nasa.gov. In this second phase, participants are requested to evaluate a surface crack test 
according to the published standard. This round robin exercise is conducted under the auspices of ASTM 
Inter-laboratory Study 732 (ILS-732). 


Objectives: 1) Determine the consistency in the interpretation of the test evaluation requirements in 
E2899. 2) Provide additional information on the analytical consistency of finite element (FE) methods as 
prescribed in the standard for future revision of the precision and bias statements. An evaluation of 
interpolated solutions as an alternative to FE will also be requested through use of the recently developed 
TASC software (Tool for Analysis of Surface Cracks). 


Requested Analysis: Interpret the provided experimental data from the 
surface crack test shown in Figure 1 in accordance with ASTM E2899, 
including the test result in K or J at the critical angle ¢,, constraint 
quantification in terms of T-stress, and interpretation of selected validity 
criteria. Spreadsheets for recording analysis results are provided to 
enhance consistency in reporting and to clarify what data are being 
requested. 


Provided Data: A compressed archive file containing the following data 
will be provided to each participant: instructions, material identification 
and stress strain data, force versus CMOD record, and the fracture surface 
photo (Figure 2) with basic crack size measurements. These data are 
intended to represent the fundamental information that would be 
provided by a mechanical test laboratory prior to evaluating the test 
result. 


Figure 1. Experimental test. 


To Participate or Ask Questions: 
Please email us: 


Douglas.N.Wells@ nasa.gov 
Phillip.A.Allen@nasa.gov 


Figure 2. Specimen Fracture Surface - All required dimensions 
are provided as an annotated figure in the electronic data 
package. See Figure 3 for complete specimen geometry. 


25 


Requested Schedule: We request that results of this Phase || assessment be returned by email (see 
“Questions?” box on other side) by February 15, 2014. As always, earlier would be very helpful. This will 
allow tabulation of the results for review at the May 2014 E08 committee meetings. 


THANK YOU for participating in this very helpful study! Doug Wells and Phillip Allen, NASA MSFC. 
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Figure 3. Specimen dimensions in inches (mm). 
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E08.07.03 Analytical Round Robin, Phase II — Precrack Measurements 
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APPENDIX B—ESTIMATION OF THE CRACK EXTENSION INITIATION ANGLE 
FOR SURFACE CRACKS WHEN NOT DISCERNIBLE 
FROM THE FRACTURE SURFACE 


B.1 Technical Background 


With the exception of E2899!, all of the current ASTM fracture testing standards such as 
E399!9 and E18202° take the 3D reality of the fracture toughness test and, through the use of vari- 
ous assumptions, simplify the problem to a 2D planar form to report an average fracture toughness 
representing the entire crack front. For surface crack fracture, the crack driving force is a strong 
function of the local crack front location, which prevents simplifying the problem into a single 
2D value. Instead, the result of the surface crack fracture test is strongly dependent on being able 
to identify the local crack front parametric angle, @. E2899 relies on post-test visual examination 
of the fracture surface to reveal the location of the initiation of crack tearing along the perimeter 
of the crack. This is generally successful for most cases when stable ductile crack extension takes 
place. The location of maximum tearing is assumed to be the site where crack extension began 
and is identified by the parametric angle, ¢;. Figure 16 illustrates this method as described in more 
detail in Annex A5 of E2899. 


Inscribed Semi-Circle of Radius a 


Figure 16. Definition of @; and the characterization of location of maximum crack tearing. 


In some cases, the visual fractographic evidence is limited with regard to identifying the 
initiation location. In the case of immediate unstable ductile fracture, i.e., fracture that occurs in 
a ‘brittle’ fashion due to a flat R-curve shape supporting no stable crack tearing, the fracture sur- 
face may not readily reveal any clear evidence of the initiation site, short of an advanced fracto- 
graphic investigation, which is outside the scope of E2899. There are also occasional cases where 
stable ductile tearing occurs in a sufficiently uniform manner around a large portion of the crack 
perimeter that a clear selection of the initiation angle is not feasible. In these cases, it is desirable 
to have a standardized method to establish the initiation angle, @,, to allow a test result to be gener- 
ated; otherwise, the test would be considered invalid and a useful result would be discarded. 
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For fracture toughness evaluations in the elastic-plastic fracture mechanics regime, the 
J-integral is evaluated using elastic-plastic finite element methods to incorporate the shifts in the 
J distribution due to plasticity. However, in most cases, initiation under uniform tension stress 
does not occur at the maximum value of the J-integral along the crack perimeter, J,, location. This 
strongly suggests that the crack front constraint conditions are influencing the angle of initiation. 
There have been many previous efforts to develop predictive methods for surface crack initiation 
using a crack front field amplitude parameter (K or J) combined in some fashion with a constraint 
parameter. Of particular interest to us is the work of Newman, Reuter, and others as described in 
references 15 and 16. The approach used in these studies was to apply a multiplicative factor to the 
local crack front field amplitude based on the local prevailing constraint conditions. Newman et al. 
characterized the local constraint using o,, a constraint parameter based on the average opening 
mode stress over the plastic zone along a ray normal to the local crack front, normalized by the 
yield strength, o,,.!>-!7.!8 The approach has been to evaluate the parameter K+q, or J+@, along 
the surface crack perimeter, with the location of maximum value of the product indicating the 
angle of initiation. The latter method shows significant improvements to the prediction of @; 
compared to using only the location of Jy: 


Based on these observations, for cases when the fracture surface cannot readily reveal the 
location of initial crack extension, a decision was made to incorporate some of this logic into 
E2899 to improve the estimation of @,. For simplicity purposes, the test standard uses the linear- 
elastic T-stress as the descriptor of crack front constraint as a function of 9, Ty. Values of Ty are 
tabulated in the standard for the requisite range of crack shapes and sizes. Given that this informa- 
tion is currently available and utilized in the standard, the motivation is to employ 7; in a similar 
predictive framework as discussed above, 1.e., providing a multiplicative influence factor for J such 
that their product estimates the crack perimeter location where driving force and constraint com- 
bine to provide the highest likelihood for crack extension. 


The ratio of T/o,,, is a good first-order estimate of the local crack front constraint condi- 
tion and is capable of predicting the constraint-induced variations in crack tip opening mode stress 
(6,,) for surface cracks under elastic-plastic conditions.*!7? Though o,, is likely not the sole influ- 
ence factor for constraint effects in ductile fracture, the usefulness of o,,, demonstrated through oy, 
indicates that 7/ ys should be capable of providing a suitable influence factor. The ratio 7/ 0,, can 
easily be related to Oy, through correlation with the parameter Q, which represents the difference 
between the constraint-influenced opening mode stress and a reference stress field, typically defined 
at T=0, evaluated at a fixed distance in front of the crack tip at ro, /J=2, where r is the radial 
distance from the crack tip.?? Figure 17 illustrates the calculation of Q as a function of T/ O,,, over 
a range of strain-hardening coefficients, m for a material with an elastic modulus, E, to Oy ratio of 
400. Note that Q(@) is defined as a difference in G5, for the surface crack compared to Os in the 
T=0 reference field, therefore, Q0=0 when T=0. To adopt the value of Q as an influence factor and 
useful multiplicative factor for J, the value of Q+1 is chosen such that, at the neutral T7=0 con- 
straint state, the multiplicative factor on Jis unity. Lower constraint values (7/ Gus 0) provide an 
influence factor <1, and positive constraint conditions present an influence factor >1. From fig- 
ure 17 it is clear that, for 7’<0, the relation 7/o,,,+ 1 is an acceptable approximation of Q+ 1. For 


T>0, T/4o,,+ 1 provides a representation of the more subtle effects of positive T-stress on o,,,. 
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Common Range of T/oy, 
for Surface Cracks 
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Figure 17. Relationship between T/o,,, and Q. 


The bilinear equation (6), provided in Annex A5 of E2899, represents a function of @ whose 
maximum identifies the estimated value for @;: 


J 
fa ees for — <0 
Pp ys ys 
Jy( 7 T 
f@)=—* Go t!| for —>0. (6) 
Dp OVs O75 


Given that 7/o,,, is a linear-elastic concept, its values are confined to the range —1< 7/o,,, 
therefore, the range of the multiplicative influence factor is from 0 to 1.25, with 1 a ae 
to the 7=0 case and having neutral influence. 


In equation (6), J, is normalized by J oe Clearly, including J, does not affect finding the 
maximum of /(@) but this normalization renders f(@) dimensionless and helps minimize the likeli- 


hood of it being misconstrued as a value with physical meaning. 


Figure 9 illustrates the a @ estimation scheme based on the J, location along with the 
currently proposed approach of finding @; by maximizing equation (6). 
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B.2 Experimental Evaluations 


A common situation where experimental results indict evaluation at @; based on the Jy loca- 
tion occur for surface cracks 1n tension with half crack surface length, c, greater than the surface 
crack depth, a, (c>a) when the J é location is often at ¢=90°. In these cases, surface crack extension 
typically initiates closer to 15°<@<35°. A small series of 2219-T8 surface crack tests are illustrated 
in figure 18. All tests are tension loading except specimen 2005-B1, which was tested in four-point 
bending. The dashed line represents perfect agreement between measured and predicted @,. Equa- 
tion (6) (open symbols) improves upon the prediction based on J f (filled symbols). 


© 155-1 
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A 155-22 
V 155-5 
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Figure 18. Performance of @; prediction based on location 
of J p Versus equation (6). 


To further investigate the efficacy of equation (6), a number of the surface crack test data 
from reference 16 were reevaluated based on available data. The test specimens are D6AC steel with 
high strength and very little strain hardening. Figure 19 shows the results of the bend tests from the 
series. The use of equation (6) provides a modest improvement for a number of the tests. Others are 
just as closely predicted by the J, location. A larger collection of reliable bend tests with known @; 
is needed. After further study, it may prove sufficient to use the J, location to estimate @; for bend 
tests in lieu of the effort required to find the maximum of equation (6). Currently, E2899 requires 
using equation (6) for all tests that do not reveal @; in the fracture surface morphology. 
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Measured Max. Tearing Location ¢ (deg) 


Filled Symbols = Jp Prediction 
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Figure 19. Performance of @; prediction for D6AC bend tests. 


Figures 20 and 21 show the results for the estimation of @; in the tension tests from the 
D6AC study. In this case, the results for J - and equation (6) are shown in separate figures for bet- 
ter clarity. Again, the most striking improvement comes in the prediction of tension tests where 
J, occurs at the crack depth. The constraint distribution moves the location of the initial crack 
extension toward 15°< @< 35°, and equation (6) generally predicts this well. 
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Figure 20. Predicted @; based on J, for D6AC tension specimens. 


Measured Max. Tearing Location ¢ (deg) 
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Figure 21. Predicted @; based on equation (6) for D6AC tension specimens. 


Figure 21 also illustrates that this predictive methodology remains an estimation that cannot 
account for all possibilities (CT-02) and there are clearly some weaknesses that may be addressed 
with further study: 


¢ Different materials exhibit a range of sensitivity to constraint due to the particulars of their 
micromechanical fracture process. The influence coefficient based on opening mode stress should 
likely vary with the constraint sensitivity of the material. 


¢ The opening mode stress model may not always be the best parameter to base the influence coef- 
ficient on. Other parameters, such as the Rice and Tracey ductile void growth model?* correlated 
to T/ O,,, May prove beneficial for some materials. 


¢ Typical material variability, local material defects near the crack front, as well as toughness 
anisotropy can strongly affect ¢;. These effects are absent from this methodology. 


* For tests with high deformation, approaching the end of the defined elastic-plastic regime in the 
standard, the linear-elastic T-stress concept becomes an increasing suspect and may not accu- 
rately reflect the true magnitude and distribution of constraint conditions around the crack 
perimeter. 


Frequently, the strong variability inherent to determining @,, whether directly from interpreting the 
fracture surface (which is often not as easy as one would hope) or from estimating ¢; using methods 
such as equation (6), is not nearly as prevalent in the calculated value of the test result (J g)> because 
the most common cause for ambiguity in @, is a nominally uniform distribution of J and constraint 
around the crack. In these cases, the J 6 result is fairly insensitive to @,. 


B.3 Conclusion 


The proposed use of equation (6) is only for cases where the fracture surface refuses to yield 
the value of @,. In these cases, rather than reject an expensive test result, it is worthwhile to have 
areasonable engineering estimate of @, for evaluating the test. The preponderance of experimental 
evidence suggests that using the @¢ location corresponding to J . is frequently inaccurate, particularly 
for tests in uniaxial tension. The previous work of Newman et al. provides a substantial foundation 
on which to reformulate the model based on the 7/o,,, constraint parameter already available as 
part of the test standard. Though not a universal panacea, the use of equation (6) generally 
represents a significant improvement over just using J, in cases where ¢; must be estimated. 
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